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Alan Coley§ and Martin Goliath| 

Abstract. Self-similar, spherically symmetric cosmological models with a per- 
fect fluid and a scalar field with an exponential potential are investigated. New 
variables are defined which lead to a compact state space, and dynamical systems 
methods are utilised to analyse the models. Due to the existence of monotone 
functions global dynamical results can be deduced. In particular, all of the 
future and past attractors for these models are obtained and the global results 
are discussed. The essential physical results are that initially expanding models 
always evolve away from a massless scalar field model with an initial singularity 
and, depending on the parameters of the models, either recollapse to a second 
singularity or expand forever towards a flat power-law inflationary model. The 
special cases in which there is no barotropic fluid and in which the scalar field is 
massless are considered in more detail in order to illustrate the asymptotic results. 
Some phase portraits are presented and the intermediate dynamics and hence the 
physical properties of the models are discussed. 



PACS numbers: 0420, 0420J, 0440N, 9530S, 9880H 

1. Introduction 

Scalar field cosmology is of importance in the study of the early Universe and 
particularly in the investigation of inflation (during which the universe undergoes 
a period of accelerated expansion [[j] ||). One particular class of inflationary 
cosmological models are those with a scalar field and an exponential potential of the 
form V{4>) = Vbe~ K ^, where k and Vb are constants. Models with an exponential scalar 
field potential arise naturally in alternative theories of gravity such as, for example, 
scalar-tensor theories, and are currently of particular interest since such theories occur 
as the low-energy limit in supergravity theories ]3|, ^| . 

A number of authors have studied scalar field cosmological models with an 
exponential potential within general relativity. Homogeneous and isotropic Fried- 
mann-Robertson- Walker (FRW) models were first studied by Halliwell || using phase- 
plane methods. Homogeneous but anisotropic models of Bianchi types I and III (and 
Kantowski-Sachs models) were studied by Burd and Barrow pj in which they found 
exact solutions and discussed their stability. Bianchi models of types III and VI were 
studied by Feinstein and Ibahez J?]], in which exact solutions were found. An analysis 
of Bianchi models, including standard matter satisfying various energy conditions, was 
completed by Kitada and Maeda || 0]. They found that the well-known power-law 
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inflationary solution is an attractor for all initially expanding Bianchi models (except 
a subclass of the Bianchi type IX models which will recollapse) when k 2 < 2. 

The governing differential equations in spatially homogeneous Bianchi cosmologies 
containing a scalar held with an exponential potential reduce to a dynamical system 
when appropriate normalised variables are defined; this dynamical system was studied 
in detail in |u|. In a follow-up paper |ll| the isotropisation of the Bianchi VII^ 
cosmological models possessing a scalar field with an exponential potential was 
further investigated; in the case k 2 > 2, it was shown that there is an open set 
of initial conditions in the set of anisotropic Bianchi VII/j initial data such that 
the corresponding cosmological models isotropise asymptotically. Hence, spatially 
homogeneous scalar field cosmological models having an exponential potential with 
k 2 > 2 can isotropise to the future. The Bianchi type IX models have also been 
studied in more detail 

Recently, cosmological models which contain both a perfect fluid and a scalar 
field with an exponential potential have come under heavy analysis [jL3| , One of the 
exact solutions found for these models has the property that the energy density due to 
the scalar field is proportional to the energy density of the perfect fluid, hence these 
models have been labelled scaling cosmologies |Q. These scaling solutions, which 
are spatially flat isotropic models, are of particular physical interest. For example, in 
these models a significant fraction of the current energy density of the Universe may 
be contained in the scalar field whose dynamical effects mimic cold dark matter. In 
Jl5[ the stability of these cosmological scaling solutions within the class of spatially 
homogeneous cosmological models with a perfect fluid subject to the equation of state 
p = (7 — (where 7 is a constant satisfying < 7 < 2) and a scalar field with an 
exponential potential was studied. It was found that when 7 > 2/3, and particularly 
for realistic matter with 7 > 1, the scaling solutions are unstable; essentially they are 
unstable to curvature perturbations, although they are stable to shear perturbations. 
Curvature scaling solutions |L6| and anisotropic scaling solutions JlTj are also possible. 
In particular, in Jl6| homogeneous and isotropic spacetimes with non-zero spatial 
curvature were studied. 

Clearly it is of interest to study more general cosmological models, and in this 
paper we shall comprehensively study the qualitative properties of the class of self- 
similar spherically symmetric models with a barotropic fluid and a non-interacting 
scalar field with an exponential potential. Self-similar spherically symmetric perfect 
fluid models with a linear equation of state have been much studied in general relativity 
[^8| ^9| . Carr & Coley ||(| have presented a complete asymptotic classification of 
such solutions and, by reformulating the field equations for these models, Goliath et 
al |2|, |22| have obtained a compact three-dimensional state space representation of 
the solutions which leads to another complete picture of the solution space. Recently, 
these models have been further studied in a combined approach P3[. The present 
analysis can be thought of as a natural extension of this recent work ]2i|, p2| , [23| . 

The Kantowski-Sachs models appear as a limiting case of the class of spherically 
symmetric models under investigation, and hence this analysis complements recent 
analyses of spatially homogeneous Bianchi models Jig] . Models with positive spatial 
curvature have attracted less attention than Bianchi models with zero spatial curvature 
or negative spatial curvature since they are more complicated mathematically. 
However, the properties of positive-curvature FRW models [l6|, Q and Kantowski- 
Sachs models || have been studied previously. 

In the next section we shall describe the governing equations of the class of models 
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under investigation. In section ^, compact variables are defined and the resulting 
dynamical system is derived in the case of spatial self-similarity. A monotone function 
is obtained. The equilibrium points and their local stability is discussed in section 
|]. The timelike self-similar case is then considered in sections |B| and ^|. The special 
equilibrium points for values of "extreme tilt" are discussed separately in section ^. 
The global results and a discussion is given in section ^. Applications in the absence 
of a barotropic fluid and in the further subcase of a massless scalar field are discussed, 
respectively, in sections ^ and [To], partially to illustrate the early-time and late-time 
behaviour of the models. 

2. Governing equations 

We shall consider spherically symmetric similarity solutions in which the source for 
the gravitational field is a perfect fluid and a non-interacting scalar field with an 
exponential potential in which the total energy-momentum tensor is given by: 

Tab = (Tpfab + T s f a b), (1) 

where 

T p {ab = ^u a u b + p (u a u b + g ab ) ■ (2) 
The perfect fluid obeys the equation of state 

P = (7-l)/*. (3) 
where 7 is a constant satisfying 1 < 7 < 2. The scalar-field contribution is given by 

1 



T B { a b 

where 



2 4>A- C + v(4>) 



9ab, (4) 



V{4>) = V Q e-*. (5) 

Since the fluid and scalar field are non-interacting, we have the following separate 
conservation laws: 

V„T p f = = V a T s f. (6) 

The spacetime is self-similar and consequently admits a homothetic vector r) a . This 
implies that the matter fields must be of a particular form. Thus, a barotropic fluid 
must have an equation of state of the form (||) |2f| . The energy-momentum tensor of 
the scalar field must satisfy: 

£*r)T s fab = Q- (7) 

This implies that 

<t> = m + \ (8) 

K 

V = e- 2 "V($), (9) 
V($) = V e- K ^, (10) 

where r\ is the variable defined by the homothetic vector ?y a , and £ is the similarity 
variable. When the similarity variable is timelike, we will use the notation t = £, x = 77. 
The homothetic vector then is spacelike, and we denote this as the spatially self-similar 
case. When the homothetic vector is timelike, we have the timelike self-similar case, 
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for which t = 77, x = £. A dot denotes differentiation with respect to the similarity 
variable. Finally, we shall define the new variable: 
1 



X = 

and for convenience we introduce the new constant 

K 



(11) 

(12) 



3. Spatially self-similar case 

In the spatially self-similar case, the line element can be written 
ds 2 = e 2x ds 2 , 

ds 2 = -dt 2 + D 2 dx 2 + Djdfl 2 , 

The kinematic quantities of the congruence normal to the symmetry surfaces are 
related to the Misner variables ((3°, (3 + ) as follows: 

9 = 3/3°, cr+ = 3/3+. (16) 

Following pljl, we will work with boosted kinematic quantities (9, a+) 



(13) 
(14) 
(15) 



(17) 



The reason for this is that it simplifies the constraint obtained from the non- vanishing 
off-diagonal component of the field equations. The metric functions B\ = D^ 1 and 
B 2 = D^ 1 then have the following evolution equations: 

B\ = ^f+Bi, (18) 
B 2 = --L(§+a+)B 2 . (19) 



The physical quantities associated with the perfect fluid are as follows: 



vT 



V v 



2/3+ -$° 



0,0 



1-v 2 



M = 



(20) 



(21) 



1 + (7 - l)v 2 

where v is the tilt variable, and fi n is the energy density along the normal congruence 
From the conservation equations (o) for the perfect fluid, we have: 



Mr, 



7 



VS[l + h-l)v 2 ] 
1-v 2 



29 



(1 - v 2 )a + + 2\/ 1 ivB 1 



• -^|l-(7-lM {7[2(7 " 1) ' + Tff+1 " 
W3 [ft - 1)(3 7 - 2)v 2 - (2 - 7 )] S t } , 

and the conservation equation for the scalar field yields 



X 



1^(26 + -„ + )X + 2kBl + V -. 



(22) 

(23) 
(24) 
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The Einstein field equations then yield the following: 
The Friedmann equation: 

H„ = - [9 2 - a\ - 3(1 + k 2 )B 2 + 3Bf - 3X 2 - 3V] . 
Constraint equation: 

o = jvfi n _ A [i + ( 7 _ i) w 2] ^ + + Sl 

Evolution equations for 6 and <7+: 

6 = - -y= [6 2 + a\ + 65+ - 3(1 + k 2 )B\ + 3X 2 - 3V 

3( 7 -l)(l-^ 2 ) 
+ l + ( 7 -l)« 2 M ' 

' ct 2 + 20ct + + 6A: 2 B 2 + 3V 



(25) 
(26) 



V3 

3 (2 - 7 ) + (37 - 2)v 2 
"2 l + (7-l)w 2 



(27) 



(28) 



From equation (25), by demanding /i„ > 0, a dominant quantity 

f = ^9 2 + 3B 2 
is identified. Thus, we define bounded variables according to 

Qa -y, Q+- -y-, Ci - ~^y—, U - —y—, W - — 

Defining an appropriate density parameter with respect to fi n , the Friedmann equation 
takes the form: 



(29) 



(30) 



3/x r , 
Y 2 



= 1 - Q% - (1 + fc 2 )C 2 - C/ 2 - VK 2 , 

while the constraint equation becomes 

= lv n n - 2 [1 + ( 7 - l)v 2 ] (Q+ + kU)C x , 

By defining a new independent variable 

, d -s/3 d 
~ dr' ~ ~Fdi' 
the evolution equation for Y 

m , 7 



F' = - <Q+ + Q 



2{Q\ + [7 2 



S2j ; 



l + (7-l)w 2 

decouples, and we are left with a reduced set of evolution equations 



Q'o = - (l - Ql) 



1-2(1-Q+-C/ 2 



7 



1 + (7 - l)w 2 



y 



Q'+ — QoQ^ 



-2(1 -Ql- U 2 ) 



7 



2fc 2 c 2 + ^ 2 + 7 1(2 ^ ) + ^- 2)w2 



1 + (7 - l)v 



l + ( 7 -l)i; 2 



(31) 
(32) 

(33) 
(34) 

(35) 

(36) 

(37) 
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C[ =2Ci 



Qo(Q 



U 2 



l-v z 



2 1 + (7- l)v 2 
{7 [2(7-l)g + 7Q+] «+ 



7 [1 - (7 " l> 2 ] 
+ [(7-l)(3 7 -2) W 2 -(2-7)]C' 1 }, 



U' 



-2(1 -Q 2 + -U 2 ) 



7 



l + ( 7 -l)v 2 



+ 2kC 2 



W 



w 2 



W = — \kQ + -U + kQ 



2(01 



IF 



7 



1 + (7 - 1> : 



(38) 

(39) 

(40) 
.(41) 



This system is invariant under the transformation 

(r, Qo, Q+, Ci,v, U, W) (-r, -Q 0) -Q+, -w, -tf , WO- (42) 

Furthermore, W > 0, and by noting the invariance under the transformation (Ci, t>) — > 
(— Ci, — u), we can without loss of generality restrict the analysis to C\ > 0. There is 
also an auxiliary evolution equation for Q n : 



a, 



7 

1 + (7 - i)u 2 



[2Q + (1 - t; 2 )Q+ + 2vd] + 2 } (43) 



In appendix Appendix A, expressions for some important fluid quantities are given. 



3.1. Invariant submanifolds 

A number of invariant submanifolds can be identified: 

• Plane symmetric: Qo = ±lj which implies B2 = (sec pH). 

• Non-self-similar Kantowski-Sachs: C\ = 0, t? = (sec |2l|). 

• Massless scalar field: W = 0. 

• No perfect fluid: fi„ = (and v decouples). 

• No scalar field: U = 0, W = and k = 0. 



The case with no scalar field was studied in 21 , where the global dynamics 
of these models was investigated in detail. The resulting state space in this case is 
effectively three-dimensional and is illustrated in that reference. We note that the 
only self-similar FRW models are the zero-curvature models which occur both as an 
equilibrium point in the plane symmetric invariant set (with Qo = ±1 = — 2Q+), and 
as an orbit in the interior of the state space. 



3.2. Monotone function 

The evolution equation for 3 = Q + + kU is given by 

,2 ^ , 



3' = Q 5 



-2(1 - Qi 



1 + (7- l)v 2 



l(2- 7 ) + (3 7 -2y 
2 l + ( 7 -l)t; 2 



a. 



(44) 
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Table 1. Summary of possible attractors for the spatially self-similar case. 



K-rings ±K 

Scalar- field dominated ±<I> 



subset of + K (-K) is always sources (sinks) 

+ $ (_$) sink (source) when 7 < §, k 2 > 1 (k 2 < 2) 



which we note is of the same form as the evolution equation for Q + in the perfect 
fluid case. Noting the form of the constraint ([33j), we consider the following function 
M (cf. 0): 

M = (1 - Qg) 2 ^S 3 ^ 4 CrV(l - « 2 )-( 2 ^. (45) 
A direct calculation then yields: 

-(37-2X2-7), 



M' 



7 v 



-(l-^)d 



(46) 



that is, M is monotonic in both the v < and the i> > regions. 

When v = the constraint yields 3Ci = 0. But since Ci = 0, i> = is an 
invariant (boundary) set, if C\ = then C\ = 0, w = always. Hence on the surface 
v = and in the interior of the state space Ci > 0. Setting v = in the evolution 
equation for v then gives 



7 



-Ci, 



(47) 



which is strictly negative. That is, all orbits in the interior region pass from v > 
to v < across the surface v = (i.e. the surface w = acts as a membrane). 
Consequently there can be no closed or recurrent orbits in the interior of the state 
space. 



4. Equilibrium points for the spatially self-similar case 

We shall display all of the equilibrium points below along with their eigenvalues. We 
will not present the corresponding eigenvectors explicitly. In what follows, e = ±1 is 
the sign of (Qo, which indicates whether the corresponding solution is expanding (+) 
or contracting (— ). The quantity Q<£ = U 2 + W 2 is the scalar-field contribution to the 
density parameter Sl n . The order of the dependent variables is (Qo, Q+, C\ 7 v,U, W) . 
The '±' suffices on the labels for equilibrium points correspond to the sign of Qo (i-e. 
the value of e) . Equilibrium points that act as attractors are listed in table El 



4.1. No scalar field (U = 0, W = 0) 
4- 1.1. K-points 

These are special points on the K-rings, defined in subsection 4.2.1. They all have 
Qo = il) Q+ = ±1: an d all other variables equal to zero. 



4-1-2. Flat Friedmann 
±F: (e,-ie,0, 0,0,0). 

n„ = f , =_o, <z pf = 2^2. 

Eigenvalues {C\ eliminated): 

1 , v 3 , , 1 , . 3 , v 37 

-(3 7 -2)e, --(2-7)6, ^(3 7 -2)e, -±e. 
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These points are saddles. 



4-1-3. Self-similar Kantowski-Sachs 

The state space contains the non-self-similar Kantowski-Sachs solutions as a boun- 
dary submanifold. In this submanifold, the self-similar Kantowski-Sachs solution 
appears as an equilibrium point. 

± KS: (f5 e '-fS f '°' ' '°)^ <2 / 3 ' 

As this solution is physical only when 7 < 2/3, we will not consider it further. 

4.2. Massless scalar field (U ^ 0, W = 0) 

4-2.1. K-rings 

±K: (e,±Vl - */ 2 , 0,0, ?7,o). 

n n = 0, n<f, = U 2 , q p{ = 2. 
Eigenvalues (Ci eliminated): 

2e, (2 - 7) [Q+ + 2e] , 2( 7 -l) e + 7 Q + , 0, 2e + Q + -p 

Each K-ring corresponds to a one-parameter family of equilibrium points (and hence 
gives rise to a zero eigenvalue). They are analogues of the Kasner solutions in the case 
with no scalar field pi] . For each K-ring, there is a subset of future or past attractors. 
+ K: sources and saddles. 
_K: sinks and saddles. 



4.2.2. M-points 

±W: (e, -k 2 fe, /, ve, kfe, 0) , / = 1/(1 + fc 2 ), 7 > §, 

- _ -7[(7-l)-^fc 2 ]±V(7-l)(2-7)(37-2)+7 2 [(7-l)-^fc^f 
Vl > 2 ~ ' (7-l)(37-2) 

Q n = 0, 00 = (iq^r) ■ 
Eigenvalues {C\ eliminated): 

-2/e, -2(l-fc 2 )/e, -(1 - k 2 )fe, F 1 ( 1 ,k)fe, 

- — [2 - 7 + (3 7 - 2)v 2 + 2 1V \ fe, 

These equilibrium points are related to the Milne points M in pi] ]. They are only 
physical (\v\ < 1) for certain ranges of 7 and k. For instance, when k < 1 it follows 
that \v2 1 > 1. When fc > 1, these points are saddles. Furthermore, examining the 
eigenvalues numerically for < k < 1, it turns out that the points always are saddles. 



4-2.3. Curvature- scaling solutions 

±X«: (2kge,-kge g,ve,ge,0), g = l/(y^yT + P), k < 1 (k 2 > 2), 

~ _ -(37-4)7-t±y/(7-l)(2-7)(37-2) + (37-4) 2 7 2 ^ 
V h2 ~ _ ( 7 -l)(3 7 -2) 
0„ = 0, =j2( 1 + fe S) ■ 

Eigenvalues (Ci eliminated): 
1 - k 2 

r — ffe, ^2(7. k,v)e, F 3 (*y,k,v)e, 

F 4 (j,k,v)e, F 5 (~f,k,v)e. 
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For k = 1 (k 2 = 2) these points coincide with ±M V , and for k > 1 (k 2 < 2) they are 
unphysical. Note that < V\ < 1 and — 1 < u 2 < only when 7 > 4/3, assuming 
fc < 1. For 7 < 4/3, it follows that \v\ > 1, and the equilibrium points are unphysical. 
Numerical evaluation of the eigenvalues shows that these points are saddles. 

4-2.4- Equilibrium lines with arbitrary v 

±$»: (e, -2^e, 0, v, ±±y/(2 - 7 )(3 7 - 2), 0) . 

fi„ = 0, ^=_(2- 7 )(3 7 - 2)/ 7 2 . 
Eigenvalues (Q+ eliminated): 

2 2 U 

0, 0, 2e, -(2- 7 )e, -e- T - 

7 7 K 

There are two zero eigenvalues for these points. The first zero eigenvalue corresponds 
to the fact that we have a line of equilibrium points. The second zero eigenvalue 
indicates that the equilibria are non-hyperbolic. For v = 0, ±1, these equilibrium 



lines coincide with the various K-rings, see subsections 4.2.1 and 7.2.1 , and these 
exceptional points mark where the K-rings change stability. The higher-order zero 
eigenvalue of ±$ v corresponds to the eigenvalue associated with the fact that ±K 
is a line of equilibrium points (and not to the eigenvalue that becomes zero due to 
the stability change of ±K), and the corresponding eigenvector is v = -^£q + + 
Perturbing the equilibrium lines ±$> v along this eigenvector, we find that 

Q' + = - 2(1 -Q 2 + -U 2 )eQ +1 (48) 
U' = - 2(1 - Q 2 + - U 2 )eU. (49) 

This is precisely the dynamical system restricted to the invariant set Qq = e, C± = 0, 
W = 0, Q n = 0. We can explicitly integrate equations ( |48"| ) and (^9|). It follows that 
Q+ is proportional to U, and the orbits in the (Q+, U) plane consist of straight lines 
through the origin with additional equilibrium points at ±& v (where Q + — — 
which are thus non-linear saddles. 

4-3. Scalar field with potential (U 7^ 0, W 7^ 0) 

There are a number of solutions with a non-zero potential listed below. There are also 
equilibrium points ±L V with variable values 

/ n 1 * n /~^\ * _ -(7-l)7V , T :r fc^^V(7-l)(2-7)(37-2) + (7-l) 2 7 2 (l~fc 2 ) 

\ 6 > U ' Vl-fc 2 ' ' ' V J ' b2 - _ ^( 7 _i)(3 7 _2) 

but these points are unphysical, since either C\ or W is imaginary. 

4-3.1. Scalar- field dominated solutions 

±$: (e,-ie,0,0,^e,^V3fc 2 -l), k 2 > | {k 2 < 6). 
Q n = 0, % = I q p[ = ^(\-k 2 ) = -1(2 - K 2 ). 
Eigenvalues (constraint degenerate): 

k 2 - 1 3fc 2 - 1 3fc 2 - 1 k 2 -1 
-e, — — e, — n — e, — e, 



k 2 ' 2fc 2 2fc 2 2fc 2 

37-4 37/fc 2 - 2 

When f2„ = 0, the constraint defines two hypersurfaces (C% = and Q + + kU = 0), 
and these hypersurfaces coincide at ±$. It turns out that the constraint is degenerate 
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(VG = 0) at these equilibrium points. Consequently, all eigenvector directions are 
physical there, so we need to keep all six eigenvalues. For k 2 — 1/3 (k 2 = 6) these 
points coincide with points in the K-rings, and for k 2 < 1/3 (k 2 > 6) they are 
unphysical. 

+ <I>: sink when 7 < |, k 2 > 1 (n 2 < 2); saddle otherwise. 

source when 7 < 4, k 2 > 1 (k 2 < 2); saddle otherwise. 
We note that corresponds to the flat FRW power-law inflationary solution p7|, 

4-3.2. Curvature-scaling solutions 

±3: (5^6, -Age, 0, 0, 5 e, ^) , g = \/(y/2y/Y+W), k > 1 (« 2 < 2). 

o — n O — 2+fc 2 _ l-fc 2 _ K 2 -2 , n 

Si„ — u, Ji0 — 2 (l + fe2) i </pf — 2 +fc2 — 2(«2 + l) ^ U - 

Eigenvalues (constraint degenerate): 



1 r_, ,01 1 



fc [2(7-1) + 7 fc 2 ] 5 e, -— [fc 2 + l± V(fc2 + l)(9A;2-7) 

-I(fc 2 + 1) 5£ , -I(fc 2 -l) ffe , i[ 7 _(2- 7 )(fe2 + l)] ffe . 

For fc = 1 (k 2 = 2) these points coincide with and for fc < 1 (k 2 > 2) they are 
unphysical. When physical, these points are always saddles. 

4-3.3. Friedmann scaling solution 

±FS: (e,-ie,0,0,f 7 e,f ^7(2-7)), k 2 < £ ( K 2 > 3 7 ). 

n n = |(2 - 3 7 _fc 2 ), = § 7 £ 2 , <? pf = 

Eigenvalues (C\ eliminated): 

i(37 "2)e, ™(2- 7 )e, 1, 



(2 - 7) ± v/(2 - 7)(127 2 F-97 + 2) 

For fc 2 = (k 2 = 37) these points coincide with ±<£, and for k 2 > ^ (k 2 < 37) they 
are unphysical. When physical, these points are always saddles. 

5. Timelike self-similar case 

In the timelike self-similar case, the line element can be written [ffsj 

ds 2 = e 2t dS 2 , (50) 

dS 2 = - D\dt 2 + dx 2 + D^dfl 2 , (51) 

D X =e fi °- 2f3 \ D 2 = e^ + . (52) 

The kinematic quantities of the congruence normal to the symmetry surfaces are 
related to /?+) by: 

9 = 3/3°, cr+ = 3/3+. (53) 

Note that for the timelike self-similar case, the symmetry surfaces x = constant are 
timelike, so the normal congruence is spacelike. As for the spatial case, it is convenient 
to boost the kinematic quantities in order to simplify the constraint |22]] : 

9 =^=(20 + a+), o-+ = ^=(e + 2a+). (54) 
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The evolution equations for the metric functions B\ = D^ 1 and B 2 = D^ 1 are 

B x = -^f+Bi, (55) 

i? 2 = --^(d + <r + )B 2 . (56) 

In terms of the coordinates we use, the physical quantities associated with the perfect 
fluid are given by 

"",«,0,0), (57) 



1 — u 2 
1 + (7 - l)w J 

where u is related to the tilt, and fit is the energy density with respect a congruence 
projected onto the surfaces of symmetry. The conservation equations for the perfect 
fluid yield 

1 - u 2 

" = V^-(7-l)] M2(7-l)g + 7* + ]« 

+^[( 7 -l)(3 7 -2)-(2- 7 ) U 2 ] J B 1 }. (59) 

There is also an evolution equation for /i t , but as it is rather lengthy and not used 
elsewhere, we will not give it here. The conservation equation for the scalar field is 

X = - -={26 + <j+)X + 2kB 2 - y , (60) 
v 3 k 

and the Einstein field equations give: 
The Friedmann equation: 

IH = V + 2 {l r l) f, - 4 - 3(1 + k 2 )B 2 - 'SB 2 - 3X 2 + 3V)(61) 



3h 2 + (t-1) 

Constraint equation: 



= ju fH - A [1 + ( 7 - l)u 2 ] (a+ + VikX)B 1 . (62) 
v3 



Evolution equations for 9 and <7+: 



■-j= (9 2 + a\ + 05+ - 3(1 + k 2 )B\ + 3X 2 + 3V 
3( 7 -l)(l- M 2 ) \ 

■ i + ( 7 _i)„ 9 **) ' ( 63 ) 

r+ +6/ 
3 (37-2) + (2- 7 )u 2 



-p (ct 2 + 26>o-+ + 6fc 2 S 2 - 3V 
v3 



'2 l + ( 7 -l). 2 (64) 
From ( |6l| ) and by demanding > it follows that 

f = \JP + 3V (65) 
is a dominant quantity. Consequently, we define bounded variables as follows: 

«» = f ft-^L, a-^s. e = ^ee) 
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The Friedmann equation becomes an equation for the density parameter fi t : 



r 2 

1 + (7 - l)w 2 



■~ w 2 + ( 7 -i) 

while the constraint becomes 

= 71A - 2 [1 + ( 7 - i)w 2 ] (g+ + fcc7)Ci. 

By defining a new independent variable 

, _ d _ V3 d 
~ d£ ~ f dx' 
the evolution equation for Y 

Y'= - {q [l + Q\ + Q Q+ - (1 + fc 2 )C? + U 2 



(67) 
(68) 

(69) 
(70) 



(7-l)(l-» 2 ) 
1 + (7 - l)u 2 



(71) 



decouples, and we obtain the following reduced set of evolution equations: 



Qo = -(1-O0) 



1 ^ ~ 



1 + Q% + QoQ+ - (1 + fc 2 )Cf + C/ 2 - ^Oot/- 



(7-l)(l-^ 2 ) 
1 + (7 - l)u 2 



(72) 



Q'+ 



Q+ j (Qo + Q+)(i - QoQ+) - ^U(i - Qo) 



+Qo 



(1 + fc 2 )C? - L> 2 



( 7 -l)(l- M 2 ) 
1 + (7 - l)u 2 



2. 2 C 2 -(l-Q 2 ) + i- (3 ^ 2) + (2 -^ 2 



1 + (7 - l)w 2 

Ci = c*i {q + Q+ + Q [q+ + QoQ+ - (l + k 2 )Cl + u 



(7-l)(l-^ 2 ) 
l + ( 7 -l)u 2 

1-m 2 



+ ^(i-QS) 



M = 



7 [u 



{ 7 [2(7-l)g + 7Q+]" 



a = 



(7-1)] 

+ [(7-l)(3 7 -2)-(2- 7 ) U 2 ](7 1 }, 

t> {-(Qo + Q+) + Qo [q+ + QoQ+ - (1 + k 2 )Cf + U 2 

(j-l)(l-u 2 ) 
1 + (7 - l)u 2 

C 2 {-Q+ + Qo 

(j-l)(l-u 2 



+ 2kC 2 --(l-Q 2 Q )(l-U 2 ), 



Q+ + QoQ+ - (1 + k 2 )C 2 + U 2 - 



l + ( 7 -l)u 2 
This system is invariant under the transformation 

(£, Qo, Q+,C u u, U, C 2 ) - (-£, -Qo, -Q+, Ci, -«, -C/, <5 2 ) 



(73) 
(74) 

(75) 
(76) 

(77) 
(78) 
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Table 2. Summary of possible attractors for the timelike self-similar case. 



K-rings ±K 

Scalar- field dominated ±$ 



subset of + K (_K) is always sources (sinks) 
+ $ (-$) source (sink) when 



Furthermore, noting the invariance under (C±, u) — > (— C\, —u) and under C2 — > — C2, 
we can without loss of generality restrict the analysis to C\ > 0, C2 > 0. 

Note that the denominator of the evolution equation for it ( |75| ) is zero when 
u = ±V7 — 1- The only way to pass this sonic hypersurface without introducing a 
shock wave is when the numerator also is zero. This defines a submanifold of the sonic 
hypersurface, and this submanifold is only physical for the u = —^/"f — 1 case. This 
severely restricts the global dynamics fl22j . 



5.1. Invariant submanifolds 

A number of invariant submanifolds can be identified: 

• Massless scalar field: V = 0, which implies Qq = ±1. 

• Static: Ci = 0,u = (see ||). 

• Plane-symmetric: C2 = (sec p2||). 

• No perfect fluid: f2 t = (and u decouples). 

• No scalar field: U = 0, Qq = ±1 and /c = 0. 

The global dynamics of the submanifold with no scalar field has been studied 
previously in p2| ], where state-space diagrams can be found. 

5.2. Monotone function 

As for the SSS case, it is possible to generalise the monotone function for the perfect 
fluid TSS case |2^| to the case of a perfect fluid with a scalar field. This is done by 
replacing S + with Q + + ktj in equation (32) of p2j . Thus, 

M = (Q+ + kUf^Cpu- 2 ^- 1 ^! - u 2 )-( 2 -^C 2 2(2 - 7) , 

with 

-(3 7 -2)(2- 7 ), 



M' = 



7« 



-(1-OC1 



M, 



(79) 
(80) 



is monotonic in the regions u > and u < 0. Furthermore, 



3 7 



7 



(81) 



which is strictly negative. Consequently, as in the SSS case, there can be no closed or 
recurrent orbits in the interior of the state space. 
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6. Equilibrium points for the timelike self-similar case 

We shall display all of the equilibrium points below along with their eigenvalues. We 
will not present the corresponding eigenvectors explicitly. In what follows, e — ±1 is 
the sign of Qo, which indicates whether the corresponding solution is expanding (+) 

or contracting (— ). The order of the dependent variables is [Qo, Q+, Ci,u, U, Ca 

The '±' suffices on the labels for equilibrium points correspond to the sign of Qo (i-e. 
the value of e) . The quantity = U 2 + Q\ — 1 indicates the presence of a non-zero 
scalar field. Equilibrium points that act as attractors are listed in table ^. 

6.1. No scalar field (U = 0, V = 0) 

The state space contains a number of solutions with no scalar field, as presented 
below. There is also a solution with variable values ^e, — gj^Ezf] e ; 0; 0; 0; 0j i $h = 

~ 16(7— l)"^ 6 ' ^ S ^°* s so ^ u ^' 10n ^ s physical only when 6/7 < 7 < 1, we will not consider 
it further. 



6.1.1. K-points 

These are special points on the K-rings, defined in section 6.2.1 
Qo = ±1, Q + = ±1, and all other variables equal to zero. 



They all have 



6.1.2. Static solution 
±T 



nanc soiuuon 

(e V^+4( 7 -i) \ 

i e ' z 3 7 -2 e ' u ' u ' u ' 37-2 y 



Eigenvalues (Ci eliminated): 



2-7 27 1 Ji 1 - 447 + 36 

e, -e, -~e± 



' 3 7 -2 ' 3 7 -2 ' 2 2(37-2) 
These points are always saddles. 



6.1.3. Regular centre 
±C°: (e,0,0,0,0,l). 

(k = 0, ft* = 0. 

Eigenvalues (constraint degenerate): 

-2e, -e, -e, e, 2e, 2e. 
These points are always saddles. 



Massless scalar field (U ^0, V = 0) 

6.2.1. K-rings 

±K: (e,±Vl-& 2 ,0,0, U,0). 

fit = 0, ft = p> 2 . 

Eigenvalues (Ci eliminated): 

e, -2e 4e + 2Q + --[/, 0, 4e+^—Q+. 

7 — 1 fc 7 — 1 
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Each K-ring corresponds to a one-parameter family of equilibrium points (and hence 
gives rise to a zero eigenvalue). They are analogues of the Kasner solutions in the case 
with no scalar field p^ j. For each K-ring, there is a subset of future or past attractors. 
+K: sources and saddles. 
_K: sinks and saddles. 



6.2.2. M-points 

±M*: (e, -k 2 fe, /, ue kfe, 0), / = 1/(1 + fc 2 ), 7 > f , 

7(Tl)- 2 f 1 7'= 2 ±\/(7-l)(2-7)(37-2)+7 2 [h-l)-V^f 
"1,2 - 2^ ■ 

fi t = o, ^ = / 2 fc 2 . 

Eigenvalues (C\ eliminated): 

-2/e, -(l-fc 2 )/e, -2(l-fc 2 )/e, F 6 ( 7 ,fc,u)e, 

[3 7 - 2 + (2 - 7 )u 2 + 2 7 ul fe. 

71/ 

These equilibrium points are related to the Milne points M in p5j. They are only 
physical < 1) for certain ranges of 7 and fc. For instance, when fc > 1 it follows 
that |it2 1 > I- When fc > 1, these points are saddles. Furthermore, for u — U2 they 
are saddles even when k < 1. Examining the eigenvalues numerically for < fc < 1, 
there are values of 7 and fc for which the u = u\ points act as attractors. However, 
this only occurs when u\ > ^7 — 1, i.e., when the corresponding equilibrium point is 
beyond the sonic hypersurface located at u = ^7 — 1. As it is impossible to cross this 
sonic hypersurface in a regular way, the +M points will not affect the dynamics of the 
models we arc interested in, even though the points are attractors for some values of 
7 and fc. This situation is also present in the case without a scalar field p2j . 



6.2.3. Curvature- scaling solutions 
±X«: (e, -±e, ^,ue, ^e, -^-Vk^l) , fc > 1 ( K 2 < 2), 

^7fc±V^- 1 )( 2 -T)(37-2) + (^) 2 7 2 fc 2 
"1,2 - — — ■ 

fit = 0, fl<f, = jp-. 
Eigenvalues (Ci eliminated): 

fc 2 - 1 



1 r 



fc ± y/l - 3fc 2 



fc 2 ' 2fc 

37-2+(2- 7 )u 2 



1 



2-yku 



For fc = 1 (k 2 = 2) these points coincide with ±M", and for fc < 1 (k 2 > 2) they are 
unphysical. Noting that fc > 1, these points are always saddles. 



6.2.4- Equilibrium lines with arbitrary u 

± $« : ( e , -2?fe, 0, u,±^(2 - 7 )(3 7 - 2), 0) . 
O t = 0, ^ = (2-7)(3 7 - 2)/ 7 2 . 
Eigenvalues (Q+ eliminated): 

2 4 2f7 

0, 0, e, -(2- 7 )e, -e-—. 

7 7 fc 
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There are two zero eigenvalues for these points. The first zero eigenvalue corresponds 
to the fact that we have a line of equilibrium points. The second zero eigenvalue 
indicates that the equilibria are non- hyperbolic. For u = 0, ±1, these equilibrium 
lines coincide with the various K-rings (see subsections 6.2.1 and 7.2.1) and these 
exceptional points mark where the K-rings change stability. The higher-order zero 
eigenvalue of ±$ u corresponds to the one that indicates that ±K is a line of equilibrium 
points (and not to the one that becomes zero due to the stability change of ±K). The 
corresponding eigenvector is v — ~^~&q + + &u- Perturbing the equilibrium lines ±<& u 
along this eigenvector, we find that 

Q' + = - 2(1 -Q\-U 2 )eQ + , (82) 
U' = - 2(1 - Q+ - U 2 )eU. (83) 

This is precisely the dynamical system restricted to the invariant set Qo = e, C\ = 0, 
C*2 = 0, fit = 0. We can explicitly integrate equations (|8^) and (|8^). It follows that 
Q+ is proportional to U, and the orbits in the (Q+, U) plane consist of straight lines 
through the origin with additional equilibrium points at ±$" (where Q+ = —2^-^-Qq), 
which arc thus non- linear saddles. 



6.3. Scalar field with potential (U ^ 0, V =/= 0) 

There is a number of solutions with non-zero potential listed below. There are 
also equilibrium points ±3 with variable values (4e, —k 2 he, 0, 0, khe, vl — A; 2 ) , h = 
1/Vl + k 2 , k 2 < 1 (k 2 > 2), but these points are unphysical since |Qo| > 1 when 
k > 0. 



6.3.1. Scalar-field dominated solutions 

±$: (2khe, -khe, 0, 0, he, 0), k 2 < § (k 2 >6),h= 1/yTTP. 
fl t =0, f2 = (2~3k 2 )h 2 . 
Eigenvalues (constraint degenerate): 

37-4, l-3fc 2 , l-3fc 2 1-fc 2 1-k 2 , 
khe, he, he, — ; he, — : he, 

7 — 1 k k k k 

2( 7 -l)- 7 fc\ 
(7-1)* 

As in the timelike case, since the constraint is degenerate we must retain all six 
eigenvalues. For k — 1/3 (k 2 — 6), these points coincide with points in the K-rings, 
and for k 2 > | (k 2 < 6) they are unphysical. 

source when 7 < 4/3 and k 2 < 2^^— =■ (n 2 > -^j); saddle otherwise. 
_$: sink when 7 < 4/3 and k 2 < 2 X ^- (n 2 > -^); saddle otherwise. 

6.3.2. Poten tial-d ominated solutions 

±Z"*: (vS^'TiTF'"*' ' )' fc< 1 ( k2 > 2 )' 

* _ (7-1)7%^^^ V(7-l)(2~7)(37-2) + (7-l) 2 7 2 (l-fc 2 ) 
U l,2 - 2^ • 

Qt = 0, = 2 ^ k i . 

Eigenvalues (C\ eliminated): 

-Qo, -(1±\/5)Q , F 8 ( 7 ,fc), 
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Table 3. Summary of possible attractors with extreme tilt. 



K-rings ±K ± 
M-points ±M ± 

H-lines ±H=f 

Scalar-field dominated ±<f> ± 



subset of + K (_K ) is always sources (sinks) 

+ M ± (_M ) sinks (sources) when 

7 >^,fc<l( K 2 >2) 

subset of +H~ (-H + ) is always sinks (sources) 

SSS: + <f> ± (- < f )± ) sinks (sources) when 

7 > f , k > 1 (k 2 < 2) 

TSS: + $ ± (-$ ± ) sources (sinks) when 

7 < I, k 2 < i (k 2 > 6) 



2(1 - u 2 ) + 7 (w 2 - 2-y/l + k 2 Q a u - 3) 
juy/1 + k 2 

For these solutions, the potential is non-zero, since |Qo| < 1- In the context of this 
paper, these solutions may be unphysical. They are on the boundary and correspond 
to non-self-similar solutions with a cosmological constant. They are always saddles. 



7. Equilibrium points at extreme tilt 

In addition to the above equilibrium points, there is a number of equilibrium points for 
which the tilt is extreme, i.e., v = ±1 or u = ±1. These are artifacts of the particular 
approach that we have adopted, and signify that the coordinates break down. These 
points are still important since orbits that are asymptotic to them may pass between 
the spacelike and the timclikc self-similar regions (at least at non- vacuum equilibrium 
points). Furthermore, for submanifolds where the tilt variable is not specified (e.g. the 
fluid vacuum submanifold) , some of these solutions are indistinguishable from similar 
ones with non-extreme tilt. Equilibrium points in this class that act as attractors are 
listed in table In what follows, the equilibrium points will be given both in the SSS 
variables (Qo, Q+, Ci,v, U, W) and in the TSS variables (Qq, Q + , Ci,u, U, 0%). 



7.1. No scalar field 



7.1.1. K-points 

These are special points on the K-rings, defined in section 7.2.1. They all have 
Qo = Qo = ±1> Q+ = Q+ = ±l,u = u = ±l, and all other variables equal to zero. 



7.1.2. C points 

In the TSS case there are equilibrium points that resemble the regular centre ±C°, 
but have u = ±1. 
±0+: (e, 0,0, ±1,0,1) (TSS). 

n t = o, = o. 

Eigenvalues (C2 eliminated): 

7-1 e±3 
£ ' ~ £ ' -^—^ — ■ 
These points are always saddles. 
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7. 2. Massless scalar field 

7.2.1. K-rings 

+ K± _K±: (e, ±Vl - U 2 , 0, ±1, U, 0) (SSS and TSS). 
O = 0, = C/ 2 . 

Eigenvalues (SSS: Ci eliminated, TSS: C\ eliminated): 

SSS:2 £ , 2(e + Q + ), _ 2 1(7 - 1)* + 7<5+ 



0, 2e + Q+ 



TSS:e, 2(e + Q+), -2 



0, 2l2 e + Q + -^ 



2- 7 
E/ 
¥' 

2( 7 -l)e + 7 Q + 



2-7 



For each K-ring, there is a subset of future or past attractors. 
+ K ± : sources and saddles. 
-K^: sinks and saddles. 

7.2.2. M-points 

a) ± M ± : (e, ~k 2 fe, /, e, kfe, 0), / = 1/(1 + k 2 ) (SSS and TSS). 



2 



O — 0, fid, — i jqjp 
Eigenvalues (SSS: Ci eliminated, TSS: Ci eliminated): 

_ 2/e , _ 4/e , _, g±^3g±jQ fc 

2-7 

-2(l-fc 2 )/e, -(l-fc 2 )/e. 
+M+: sink when 7 > ^ppp , fc < 1 (k 2 > 2); saddle otherwise. 
_M~: source when 7 > 2 g 3 _^ - , fc < 1 (k 2 > 2); saddle otherwise. 

b) ±M=F : (e,-k 2 fe, f,~e,kf e,0), f = 1/(1 + k 2 ) (SSS and TSS). 

o = o, n^=(^) 2 . 

Eigenvalues (SSS: C\ eliminated, TSS: C\ eliminated): 

-2/e, -(l-fc 2 )/e, -2(l-fc 2 )/e, -2(1 - k 2 )fe, 0. 
The zero eigenvalue is due to the fact that these points are the end points of the 



equilibrium lines ±H (see section 7.2.3). 

+ M~: sink when k < 1 (k 2 > 2); saddle otherwise. 

_M + : source when k < 1 (k 2 > 2); saddle otherwise. 

7.2.3. H-lines 

±HT : ( e , Q + , 1 + eQ+,-e, k(l + eQ+)e, 0), -1 < eQ+ < (SSS and TSS). 

O = -2e(l + eQ+) [(1 + k 2 )Q+ + ek 2 ] , % = k 2 {l + eQ+) 2 . 
Eigenvalues (SSS: C\ eliminated, TSS: C\ eliminated): 

-2(l + 2eQ+)e, -2(1 + 2eQ+)e, -(l + 2eQ+)e, 

-2(l + eQ+)e, 0. 
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These equilibria consist of lines of equilibrium points. The eigenvector direction along 
the lines is v — eq + + kejj. The end points of the lines are the M-points ^M^ at one 
end and points of stability change on the K-rings ±K T at the other end. 
+H~: always contains at least a subset of sinks. Solely sinks when k > 1 (n 2 < 2); 
sources and saddles otherwise. 

_H+: always contains at least a subset of sources. Solely sources when k > 1 (k 2 < 2); 
sources and saddles otherwise. 



7.2-4- Curvature- scaling solutions 

+ X± _X ±: (2kg e,-kge,g,±l,ge,0), 
g = l/(y/2y/T+¥), k < 1 {n 2 > 2) (SSS), 
(e, -±e, ± ±1, fe, ^ VW=l) , k > 1 (n 2 < 2) (TSS). 
Q, = 0, Sfy = 2 (i+A;2) • 

Eigenvalues (SSS: Ci eliminated, TSS: C\ eliminated): 

SSS : - 2 |^i[sgn + -^^<?e, We, 

Fi (fc)e, -Fn(fc,t;), where FiO(fc) > 0, and F n (k) > 0. 
TSS : ^(k 2 - --j-[fce+ sgn(«)], 

-I[fce+sgn( V )]fc^, -^±^^4^). 
For k = 1 (k 2 = 2) these points coincide with ±M ± . They are always saddles. 



7.3. Scalar field with potential (U ^ 0, W ^ 0) 

7.3.1. Scalar- field dominated solutions 

+ $±, _$±: ( f ,46,0,±ll t iy5M ), fc 2 > i ( K 2 <6) (SSS), 
(2fc/ie, -fc/ie, 0, ±1, he, 0), /i = 1/VT+I 2 , fc 2 < ± (k 2 > 6) (TSS). 
fi = O,fi = 3/4. 

Eigenvalues (SSS: Q+ eliminated, TSS: Q+ eliminated): 

SSS:-^e, -4^1 e> -|^ie, -2^ 2 e , -2W 2 e. 

TSS:^(l-3fc 2 ), ^(l-3fc 2 ), ^(1-fc 2 ), ^(1-fc 2 ), 

-2*p±hk. 

2- 7 

For fc 2 = 1/3 (k 2 = 6) these points coincide with points on the K-rings ±K ± . 
SSS: 

+ $ ± : sinks when 7 > |, k > 1 (k 2 < 2); saddles otherwise. 

sources when 7>|,/c>1(k 2 <2); saddles otherwise. 

TSS: 

sources when 7<|,/c 2 <l/3 (k 2 > 6); saddles otherwise. 
_$ ± : sinks when 7<|,/c 2 <l/3 (k 2 > 6); saddles otherwise. 
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7.3.2. Curvature- scaling solutions 

+5 ± ,_S ± : (^e,-kge,0,±l,ge,^), g = l/(V2VT+J?), k>l (n 2 <2) (SSS), 
(^e,-k 2 he,0,±l,khe,VT^ P), h = l/y/T+l?, k < 1 (k 2 > 2) (TSS). 

O = O , - 2+fc 2 
" — u - — 2(l+fc 2 ) ■ 

Eigenvalues (SSS: Q+ eliminated, TSS: Q+ eliminated): 

SSS : — ^ " — — _^ — e, —-rik 2 — l)qe, — — - — e, 
2V2k k K Jy 2kg 

2 [2(7-l) + (37-2)/c 2 ] 5 e. 



(2-7)* 



TSS: -ie, (l-fc 2 )/ l£ , -_L e ±IV7-9A: 2 , 

2( 7 -l)-(2- 7 )*» 

2-7 

For fc = 1 (k 2 = 2) these points coincide with ± < E >± . They are always saddles. 



7.3.3. Potential-dominated solutions 
± 



Z±: (\/i^P^ 6 ' ' TT^F' 1 * 11 ' ' )' k < 1 (k 2 > 2). 

fit = o, r2^ = 2 ^ k .± . 

Eigenvalues (C\ eliminated): 

-Qo, -d ± V5)Q , _ 2 ^±EQ^M , 

V 1 + k- 

2 (37 - 4) sgn (u) + 2( 7 - 1) yT + FQo 

(2- 7 )\/rTfc 2 

These solutions are only physical in the TSS region. They are always saddles. 
8. Global results and discussion 

Due to the existence of monotone functions and the fact that there are consequently 
no closed or periodic orbits in the physical state spaces we can obtain global results 
for the dynamics by studying the local stability of the equilibria. 

Indeed, from the monotone functions obtained in the spatially self-similar (SSS) 
case (f45[ ) and the timelike self-similar (TSS) case (|7^) we can immediately deduce from 
the monotonicity principle |26| that all orbits have Q 2 , — > 1, Q+ + kU — > or C\ — > 
(or an extreme value for v) asymptotically in the SSS case (and similarly in the TSS 
case). Moreover, we can also see immediately that by setting the right-hand-sides of 
equations (38) and (|4l| ) to zero that either C\ = or W = 0, or if both are non-zero 



then kQ+ +{7 = 0; this latter case yields very severe constraints on any possible 
equilibrium points. In fact, from the local analysis of the equilibria we can determine 
all of the sinks and sources. In both the SSS and TSS cases a set of massless scalar 
field models lying on the + K-ring act as sources (i.e., early-time attractors) and a 
set of massless scalar field models lying on the _K-ring act as sinks (i.e., late-time 
attractors), and for certain ranges of the parameters (e.g., n 2 < 2) the equilibrium 
point + $ with U ^ and V / 0, corresponding to the ever-expanding inflationary 
flat FRW power-law solution are sinks (i.e., late-time attractors). [We 
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note that the equilibrium point _<!> corresponding to e = — 1, which acts as a source, 
represents an ever-contracting solution and therefore is of less physical importance, 
although it does serve to classify all of the possible orbits in the state space.] Hence 
we have the global results that models that are initially expanding always expand 
from an initial singularity and always recollapse to a second singularity (when k 2 > 2) 
or either recollapse or expand forever towards a flat FRW power-law solution (for 
k 2 < 2). This global behaviour is the same as that for positive curvature FRW models 
and Kantowski-Sachs models |^9| and for Bianchi type IX models || [| . Models that 
expand from an initial singularity and recollapse to a second singularity are said to 
satisfy the positive-curvature recollapse property J3(], |3lJ . Models that expand towards 
the flat FRW power-law solution isotropise and inflate to the future and are said to 
satisfy the cosmic no-hair theorem g|. The time- reverse of the above dynamics is 
also possible (essentially e — * — e; although we have included this in the analysis these 
models are of less interest physically). Solutions in which the shear and the kinetic 
energy of the scalar field dominate are analogues of the Kasner and Jacobs solutions 
[p6| , and a rigorous study of the structure of the singularity, which is non-oscillatory, 
for a general class of analytic solutions of the Einstein field equations coupled to a 
massless scalar field has recently been presented [^2j . 

However, there are some aspects of the global dynamics of the self-similar, 
spherically symmetric models that are different. The complete set of attractors for 
different values of the parameters 7 and K are summarised in table ^. Some of these 
differences are quite subtle. First, we note that the flat FRW power-law inflationary 
solution corresponds both to a set of non-tilted equilibrium points ±$ and to points 
at extreme tilt i"!^. There is a bifurcation of the equation-of-state parameter at 
7 = 4/3 in that the power-law inflationary solution $ is a non-tilted attractor ±$ for 
7 < 4/3 and an extreme-tilt attractor ±$> ± when 7 > 4/3. For 7 = 4/3 there exist 
lines of equilibrium points with arbitrary tilt v (or u in the TSS case). This type 
of 7-dependent behaviour has also been found in Bianchi type V two-fluid models 
p3| |3~4j . Second, there exist additional M-point attractors (for 7 > ^rxrp , n 2 > 2) 
at extreme tilt. The significance of these is less clear, although they are important for 
the matching of orbits and they are related to critical phenomena |35| . We shall discuss 
this further in sections ^ and [h]. We recall that solutions in the SSS region and the 
TSS region can be matched across the surface of extreme tilt via the equilibrium points 
with extreme tilt (see earlier and |2l], f22|). In addition, a comprehensive analysis of 
the matching of solutions would be necessary in order to obtain a complete knowledge 
of the intermediate dynamics of the models. Clearly the intermediate behaviour of 
the models under investigation will be quite different to that of the models previously 
studied. 

We note that all of the equilibrium points with non-negligible matter, namely 
the non- vacuum Flat Friedmann ±F and the Friedmann Scaling ±FS equilibrium 
points, are saddles. This means that the perfect fluid is not dynamically important 
asymptotically. In order to understand the asymptotic behaviour of the models we 
consequently need only study the vacuum models (i.e., the invariant boundary with 
Q„ = in the SSS case and the invariant boundary with f2 t = in the TSS case). 
We shall discuss the fluid vacuum models further in section |j| The matter will play 
an important role in describing the dynamics of the models at intermediate times, 
and hence the physical properties of the models. We shall illustrate some of the 
intermediate dynamics in the next two sections. However, in the SSS case we know 
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Table 4. Summary of possible future attractors. Past attractors are obtained 
by changing the sign of Qo, which is indicated by the lower- left index. This is 
also the sign of the expansion. Comments: (1) The attractors +H _ have extreme 
tilt; orbits may pass between the SSS and TSS regions there. (2) For 7 = 4/3 
the SSS power-law solution switches from non-tilted + <3> to +"3?^ at extreme tilt; 
precisely at the bifurcation value there appear lines of equilibrium points between 
these two. (3) In the TSS case, the non-tilted _$ and the tilted coexist 
for certain parameter values. However, they are always separated by the sonic 
hypersurface. (4) In some regions further conditions must be fulfilled for some 

attractors: *) The attractor _|_M _ must satisfy 7 > \ K i+2 ' wnere ( as K ^) 

g 2(3k 2 +2) <!•**) The n on-tilted TSS attractor must satisfy k 2 > 





1 < 7 < 6/5 


6/5 < 7 < 4/3 


4/3 < 7 < 2 


K 2 < 2 


_ K-rings 
+H _ -line 

SSS power-law +<E> 


_ K-rings 
+H~-line 

SSS power-law + $ 


_ K-rings 
+H~-line 

SSS power-law + $ ± 


2 < k 2 < 4 


_ K-rings 
+ H~-line 


_ K-rings 
+ H~-line 

SSS M-point +M- * 


_ K-rings 
+ H _ -line 

SSS M-point +M" 


4 < k 2 < 6 


_ K-rings 
+ H~-line 


_ K-rings 
+ H~-line 

SSS M-point +M- * 
TSS ** 


_ K-rings 
+ H _ -line 

SSS M-point +M- 


K 2 > 6 


_ K-rings 
+ H _ -line 
TSS _$ ** 
TSS _$ ± 


_ K-rings 
+ H~-line 

SSS M-point +M- * 

TSS 

TSS 


_ K-rings 
+ H _ -line 

SSS M-point + M~ 



from the behaviour of the monotone function that the subcase C\ = is important 
asymptotically (in the TSS case the analogous case is C\ = 0, leading to the static 
models p2]| ). Moreover, when fl n = 0, the constraint ( p3[ ) leads to (Q+ + kU)C\ = 0. 
Clearly the invariant set C\ = (and v = 0), corresponding to the (non-self-similar) 
Kantowski-Sachs models, is of vital importance, and knowledge of the dynamics of 
the Kantowski-Sachs models is crucial for a complete understanding of the dynamics 
of the models under consideration here. In addition, all of the interesting transient 
dynamics with non- negligible matter (e.g., the non- vacuum Flat Friedmann and the 
Friedmann Scaling saddle equilibrium points, as well as the power-law attractors) 
occurs in the Kantowski-Sachs invariant submanifold. Consequently, we shall study 
the Kantowski-Sachs models in more detail elsewhere ]29| . 

We note that a set of massless scalar field models lying on the K-rings act as 
sources and sinks (i.e., early- and late-time attractors). It is therefore also of interest 
to study the self-similar, spherically symmetric massless scalar field models more fully. 
Indeed, such a study will also be of relevance in the study of critical phenomena (see 
section [n]). We shall return to this in future work. 
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K 
M 
X 
$ 



Table 5. Equilibrium points in the f2 n = 0, Q-\ 
subscripts on labels have been suppressed. 



-kU submanifold. Note that 



u 



w 



± 



e 

2 k 



k e 



1 + fc 
1 



c 

ySyT+F , 

2k ' 



2k 6 

l, e 

v / 2y / I+F' 



i 





fc < 1 (k 2 > 2) 
fc 2 > 1/3 (k 2 < 6) 
fc > 1 (k 2 < 2) 



9. Fluid vacuum 

When there is no barotropic fluid present (f2 n = and = 0, respectively), the 
constraint gives rise to two separate invariant submanifolds: either Q+ = — kU 
(Q+ = -kU), or else C t = (C*i = 0). The Q + = -W (Q+ = -A;Z7) submanifold 
is particularly interesting, as it contains all the sinks and sources of the more general 
models under consideration. Furthermore, this submanifold is three-dimensional, and 
hence lends itself to visual presentation. In addition, the submanifold C\ = will be 
studied in detail in [E9|. 



9.1. Spatially self-similar case 

In the SSS case, the reduced dynamical system becomes (eliminating the variable Ci): 
Q' - (1 - Ql) [1 - 2(1 + k 2 )U 2 ] , (84) 

- 2 (^- r n?) [i-(i+fc 2 )f/ 2 ] 

1 / 2fc 2 



U' 



k 



1 



1 + fc 2 



W 



l 



fc 2 ) [ 2Q U-- | fir. 



(85) 



(86) 



The equilibrium points of this system are listed in table pi They constitute a 
subset of the points listed in sections ^ and [7j In figures |l| andg, some examples of 
state-space diagrams for this model are displayed. 



9.2. Timelike self-similar case 

In the TSS case, the reduced dynamical system becomes (eliminating the variable Ci): 



Q' = -(i-Qo) 



1 



-(Qo-2kU)U 



(87) 



U' = - ^ ( 1 - Ql + 2kQ U - 



2fc 2 



1 + fc 2 



l-(l + fc 2 )/7 2 



QoU 



TTfc 2 " 1 2 ' 




c'=c 2 



(i-Qo 



2kQ U)- 



k 



Qo(i 



<-^9 / 



The equilibrium points of this system are listed in table 
subset of the points listed in sections ^ and ^ In figures || and 
state-space diagrams for this model are displayed. 



(89) 

They constitute a 
some examples of 



9.3. Discussion 

In both the SSS and the TSS cases there are always attractors dominated by the kinetic 
part of the scalar field, corresponding to the K-equilibrium points. Consequently, 
there will always be solutions that expand from a K-singularity and recollapse to a 
K-singularity. When k 2 > 2, the M-equilibrium points also are attractors. These 
correspond to dispersing solutions. Finally, there are ^-equilibrium points, which 
correspond to power-law inflationary solutions when k 2 < 2. In the SSS case they act 
as attractors when k 2 < 2, while in the TSS case they are attractors when k 2 > 6. 

To summarise, when k 2 < 2 the SSS case contains recollapsing K — > K solutions 
and power-law inflationary solutions K — > $. This is the same situation as for the 
Kantowski-Sachs models and closed Friedmann models examined in ]29[ ]. In contrast, 
when k 2 < 2 the TSS case contains only recollapsing solutions. When n 2 > 2, both 
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Figure 2. The state space for the Q n = 0, Q+ = —kU submanifold for k > 1 
(k 2 < 2). Dashed curves and white arrows are screened. The dotted curve 
between +M and +H is in the interior of the state space. The subscripts on the 
K-points refer to the sign of Qo and U, respectively. The state space is bounded 
by the following invariant submanifolds: the bottom is the massless scalar field 
submanifold (note that only part of the bottom is depicted - for a complete 
picture see Fig. hid), the half-disks are plane-symmetric submanifolds, and the 
half-cylinder is the Kantowski-Sachs submanifold (Ci = 0). Note that orbits and 
equilibrium points in the rear half have been suppressed for clarity. 



Table 6. Equilibrium points in the Qt = 0, Q_|_ = —kU submanifold. Note that 
subscripts on labels have been suppressed. 





Qo 


u 


c 2 




c 


e 





1 




K 
M 


e 
e 


l+k 2 t 








X 


e 


2k € 




k> 1 (k 2 < 2) 


$ 


2k 







k 2 < 1/3 (k 2 > 6) 


vT+fe 4 


Vi+k^ e 


z 


/i-k\ 
V i+k* e 








k < 1 (k 2 > 2) 



the SSS and the TSS cases contain recollapsing solutions and also dispersing solutions 
K — ► M. The asymptotics is thus similar to the fluid-only case (^3|, where the generic 
regular solutions either are recollapsing K — > K solutions or ever-expanding solutions 
K — ► M. Additionally, when k 2 > 6 the TSS case also contains dispersing K — > $ 
solutions that are non-inflationary. 

10. Massless scalar field 

Critical phenomena in gravitational collapse were first found by Choptuik fl36| in the 
study of a massless scalar field, and remain an active field of research (see, e.g., |$7| 



Self-similar models with a fluid and a scalar field 



20 




% K 



(a) 



(b) 




(c) 



(d) 



Figure 3. The state space for the Qt = 0, Q+ = —kU, C2 = submanifold for 
various values of k: (a) k 2 < 1/3 (re 2 > 6); (b) 1/3 < k 2 < 1 (2 < re 2 < 6); (c) 
fc = 1 (re 2 = 2) and (d) A; > 1 (re 2 < 2). The subscripts on the K-points refer to 
the sign of Qo and U , respectively. Note the line bifurcation occurring for k = 1 
(re 2 = 2). 



and references therein). The solution at the threshold of black- hole formation in 
spherically symmetric radiation fluid collapse, corresponding to a = 1/3, was studied 
by Evans & Coleman |38|| . In J|5| a new class of 'asymptotically Minkowski' self- 
similar spacetimes were presented, which were shown to be intimately related to the 
so-called critical phenomena which arise in spherically symmetric gravitational collapse 
calculations Q . 

Here, we present the governing equations for a self-similar massless scalar field in 
spherical symmetry. In this case, <j> is still of the form 

= $(£) + \/2fo7, (90) 

and the corresponding equations are formally obtained by setting V — 0. In the 
presence of a barotropic fluid we then essentially have a non-interacting two-fluid 
model |2(], 39 , in which the massless scalar field can be identified with a stiff perfect 
fluid and the two fluids are separately conserved. We can then deduce that the models 
evolve from the massless scalar field model to the single-perfect fluid model [^0[ p9[ . 
Hereafter, we shall assume that there is no barotropic fluid present, and that we are 
investigating a special case of the fluid vacuum model studied above. The massless 
scalar field equations without perfect fluid are then obtained by subsequently setting 
£1„ = and 0* = 0, respectively. This leads to the decoupling of v and C\ in the 




Figure 4. The state space for the Ot = 0, Q+ = —kU submanifold for k 2 < 1/3 
(k 2 > 6). Dashed curves and white arrows and circles are screened. The subscripts 
on the K-points refer to the sign of Qo and U, respectively. The state space is 
bounded by the following invariant submanifolds: the bottom is the massless 
scalar field submanifold (note that only part of the bottom is depicted — for a 
complete picture see Fig. j^a), the half-disks are plane-symmetric submanifolds, 
and the half-cylinder is the static submanifold (C*i = 0). Note that orbits and 
equilibrium points in the rear half have been suppressed for clarity. 



SSS case, and u and C\ in the TSS case, respectively. Furthermore, in the SSS 
case W = 0, while in the TSS case, it follows that Qo — e (= ±1)- The remaining 
dynamical systems are thus three-dimensional systems in (Qo, Q+, U) and (Q+, U, C2), 
respectively. The constraint leads to two separate regions: either C\ — or else 
Q + = —kU. Note that the latter is a special case of the models treated in section ||. 
Here, we briefly summarise the governing equations for these models and list all of the 
equilibrium points. 

10.1. Spatially self-similar case 

In the SSS case, the Fricdmann equation becomes 



= 1 - Q\ - (1 + k 2 )Cl - U 2 , (91) 
which implies that 

Z = 1 - Q\ - U 2 = (1 + k 2 )Cf. (92) 
The constraint becomes 

= (Q+ + k,U)Z, (93) 
and the reduced dynamical system is 

<3o = -(1-Qg)(l-2Z), (94) 

Q'+ = - 2Z (q Q+ + , (95) 

U' = _2z(q ^- T ^). (96) 
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Table 7. Equilibrium points in the Q n = 0, W = submanifold. Note that 
subscripts on labels have been suppressed. 





Qo 


Q+ u 




K-rings 




Ql + U 2 = 1 




M 


e, 


-k 2 fe kfe 


/ = l/(l + fc 2 ) 


X 


2kge 


-kge ge 


g = l/(V2Vl + fe 2 ) 



The constraint implies that either Z or Q + + kU must be zero. When Z — 0, both 
<5+ and U are constants (subject to 1 — Q\ — U 2 = 0), whereby the dynamical system 
reduces to a single evolution equation for Co, and consequently Qo is monotonically 
decreasing. If the constant values of Q+ and U do not satisfy Q+ + kU — 0, then 
the dynamics in the invariant set Z = does not intersect with the dynamics in the 
invariant set Q+ + kU = 0. The latter is contained within the fluid vacuum case 
studied in the previous section (see figure [l]). The equilibrium points of the system 
are listed in table 0. 



10.2. Timelike self-similar case 

In the TSS case, the Friedmann equation becomes 

= 1 - Q 2 + - (1 + k 2 )Cf -U 2 - Cf , 
which implies that 

Z = 1 - Q\ - U 2 = (1 + k 2 )C\ + C\. 
The constraint becomes 

= (Q+ + kU)(C 2 -Z), 
and the reduced dynamical system is 



Q'+ = -eQ+(2Z-Ci 
U' = -eU(2Z-Ci)- 



2k 2 



l + k 2 
2k 



(C\ 



-z), 

Z), 



(97) 
(98) 
(99) 

(100) 

(101) 
(102) 



l + k 2 

eC 2 {\ + C\-2Z). 

The constraint implies that either Z = C 2 or Q + = —kU. When Z = Cf, the 
dynamical system becomes two-dimensional. The evolution equation for C2 decouples 
and C2 and Q+ (or U) are monotonic. If the values of Q + and U do not satisfy 
<5+ + kU = 0, then the dynamics in the invariant set Z = Cf does not intersect with 
the dynamics in the invariant set Q+ + kU = 0. The latter is contained within the 
fluid vacuum case studied in the previous section, and corresponds to the semi-disks 
at Qo = ±1 in figure ^. The equilibrium points of the system are listed in table ||. 



10.3. Discussion 

The dynamics is different in the various invariant submanifolds. The Z = 
submanifold of the SSS case only contains recollapsing K — ► K solutions. The dynamics 
in the Q + — —kU submanifold of the SSS case is more complicated (see figure |l|); 
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Table 8. Equilibrium points in the Qt = 0, Qo = e submanifold. Note that 
subscripts on labels have been suppressed. 





Q+ u 


c 2 


c 





1 


K-ring 


Q 2 + + U 2 = 1 





M 


-k 2 fe kfe 


/ = l/(l + fc 2 ) 


X 


~2 e 2k 6 


±dr^ 2 - 1 



when k 2 < 2 there are only recollapsing K — > K solutions, when 2 < k 2 < 6 there are 
recollapsing solutions and dispersing K — > M solutions, and when n 2 > 6 there are also 
singularity- free bouncing M — > M solutions, analogous to the bouncing Friedmann - 
Lcmaitrc solutions. The TSS case only contains dispersing solutions. In the Z = Cf 
submanifold there are K — > C solutions, while in the Q + = —U submanifold there are 
K — ► M solutions. 



Appendix A. Fluid quantities 



For the non-tilted (v = 0) equilibrium points of the SSS case, we have given the 
deceleration parameter q p f with respect to the fluid congruence. The necessary 
expressions are summarised below. 

The expansion of the fluid congruence is given by 

Opt - V a u° f 

Ye~ x ( vv' - - - \ 

+ 2Q + Q+ + 3vC\ , (A.l) 



(A.2) 



v^vT^? ^ 1 _ v 2 
and the deceleration parameter, defined by 



3 P f 



1 



9 2 



is given by 



5 P f 



= - 1 



oh , , , v(2Q + Q + ) + C 1 (5-v 2 ) ^ 

Z{ 40 • T 12 



1 + 2^ , 2 



vv 



- (2Q + Q+ + v') {Q+ + Qo [2{Q\ + U 2 ) 



+ - 



7 



7^n 



|- {2(Qo-Q+) + 3vC 1 }vC 1 ) . (A.3) 



1 + (7 - 1)« 2 

Appendix B. Transformation between the SSS and TSS variables 



The two sets of variables basically differ only in the choice of dominant quantities 
(although one has to be careful, as the change of causality may result in sign change). 



Self-similar models with a fluid and a scalar field 



30 



(B.l) 



Defining 

= V^o - ^ 2 (B.2) 

=i— , (B.3) 

the relations between the variable sets become 

Qo = (Qo, Q+ = (Q+, ^i = C<5i, « = « _1 , tf = C#> (B.4) 



m> = cVQo-i 5 c 2 = ry^-i. (b.5) 
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